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Many-body quantum coherence is responsible for many 
astonishing properties of composite quantum systems, 
what concerns their spectral, e.g., specifically ground 
state pQ properties but also their quantum dynamics: 
think of superconductivity where coherently composed 
quantum systems function as carriers that guarantee loss- 
less charge transport [2), or of enhanced exitation trans- 
port in photosynthetic processes, which is under debate 
to be of quantum origin [3 . The inherent many-body 
character of such features makes them difficult to under- 
stand in terms of single particle observables which are 
the categories of our classically trained intuition. Con- 
sequently, we need a better understanding of the nature 
of many-body quantum coherence, and, a forteriori of 
multi-partite entanglement. 

So far, the structure of entangled states is still widely 
unexplored. Only for pure states of bipartite systems 
do we have an exhaustive understanding based on the 
Schmidt decomposition i\. But in multipartite systems, 
where there is no comparably useful tool, we have only 
very partial knowledge. We know that there exist two in- 
equivalent classes of entangled states in three qubit sys- 
tems [5], the GHZ and W states. Four qubit systems 
reveal even a number of distinct families, i.e. contin- 
uously parametrized sets of classes, of entangled states 
[BJ. Further efforts [7HH], still restrict to rather small 
and low-dimensional systems, which hitherto leaves us 
with a rather intransparent picture of the general struc- 
ture of multipartite entanglement. 

The defining property of entanglement is that it can 
not be created by local manipulation. Formally, this is 
expressed in terms of a Kraus operator K [T2 [T3]: if 
there is an operator K such that |$) = K\^}, then there 
is a generalized measurement, also referred to as POVM 
[T5] , in which one of the outcomes is associated with the 
transition of the state to |$). If the Kraus operator 
acting on an TV-partite system is a simple direct product 
of local Krauss operators Fi, i.e. if K = F± C9 . . . g) Fjv, 
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then the measurement can be implemented in terms of 
local measurements on the N subsystems. 

In bipartite systems a state |$) can be created from a 
state \*S?) by local manipulation if and only if the Schmidt 
rank, i.e. the number of finite Schmidt coefficients of 
|$), docs not exceed that of |^). Thus, starting out 
from a state with maximal Schmidt rank in a given fi- 
nite dimensional bipartite system, any other state can be 
obtained through local manipulation only. This is funda- 
mentally different in multipartite systems where there is 
in general no such so-called maximally entangled state, 
but there are distinct classes or families of multipartite 
states. Here, two states |$) and \^f) define two differ- 
ent equivalence classes if there are no local Kraus oper- 
ators Fi (i — 1,...,N) and Gi (i = l,...,N) such that 
|$) = F x <g> . . . <g> F/vl*) and |$) = G\ <g> . . . <g> G N \$). 
In tripartite qubit systems there are the two classes 
of \W) = (|211) + |121) + |112))/V3 and \GHZ) = 
(| 111) + |222))/\/2 states [5], and in fourpartite qubit 
systems there are even continuous families of inequiva- 
lently entangled states |BJ . Only if the dimension of one 
of the subsystems is as least as large as the product of the 
dimensions of all other subsystems there is a maximally 
entangled state in the above sense [13] . 



I. REPLICAS OF QUANTUM SYSTEMS 

The criterion [15] that we exploit here to identify struc- 
ture of many-body entangled states is based on permuta- 
tion symmetries of multiple replicas of a quantum state. 
'Multiple' in the present context means that rather than 
considering a quantum state |$) of an A^-partite quan- 
tum system, we consider M replicas of this state vector, 
i.e. the M-fold tensor product of the state vector |$) 
with itself. Such permutation symmetries must not be 
mistaken with symmetries under permutations of sub- 
systems which are frequently considered |16H19j . In or- 
der to avoid such confusion between the different sub- 
systems and different replicas of a quantum state, we 
will adopt the convention of to write the different 
replicas of a quantum state in different rows. All terms 
associated with a specific subsystem will the be within 
a column. That is, if we consider a bipartite quan- 
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turn state |*> = £V 



\j), then two replicas read 
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we will denote the four-component state vector as 
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or, in short hand notation, as 
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A. Permutations 

As we will discuss in the remainder of this article we 
can use permutations on several replicas of a subsystem 
to distinguish different classes of entangled states. In the 
case of two replicas, the (unique) permutation operator 
II associated with one subsystem is defined via 
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and, in the case of M replicas there are Ml inequiva- 
lent permutations 121) per subsystem that are defined 
similarly to Eq. In the case of ([I]) there are two dis- 
tinct permutations possible: one that interchanges 
and ^2) on subsystems one and another one that inter- 
changes \ji) and I j'2) on the subsystems two. Any sepa- 
rable state l^i) = \<pi) <X> ^2} is invariant under either of 
these permutations: 
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This is different for an entangled state, that we can 
write in its Schmidt decomposition as | = aAi|H) + 
\J\i\T2). Its duplicate version reads 
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This object is actually altered by a permutation on the 
first subsystem: 
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and similarly for a permutation on the second subsystem. 
Consequently, the projector P_ = (1 — II)/2 onto the 
antisymmetric component of Hi <£> Hi ( or analogously of 
Hi®Hi ) can now be used to probe the invariance of any 
given duplicate state I*) ® |\&) under II: any separable 
state is invariant under II and is thus mapped on the 
zero- vector 
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whereas any entangled state remains finite: 
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Later-on, we will see that similar distinctions also ex- 
ist for different classes of multipartite entangled states. 
However, before doing so, let us formulate our criterion 
to identify inequivalent entanglement classes. 



B. Criterion for the classification of many-body 
entanglement 

The central property of any permutation operator II 
acting on M replicas of a single-particle Hilbert space is 
that it commutes with the M-fold tensor product of any 
operator F: 
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Permutation operators are ideally suited to distinguish 
different classes of entangled states, since this relation 
naturally generalizes to M replicas of an iV-partite sys- 
tem, if we introduce introduce the operator 
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A is composed of permutations , where 'j' denotes 
the subsystem that the permutation is associated with, 
labels different permutations, and if i 1 '"f^ are complex 
numers. Any such operator A commutes with with M 
replicas of a local Kraus operator 
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With this at hand, we can rephrase the criterion for the 
classification of many-body entanglement [15] that we 
will exploit in the following: 



If an operator A as specified in Eq. ( 10 1 exists 
such that 
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then the state |$) cannot be obtained from |*) 
through local manipulation. 
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The proof of this assertion follows by contradiction: 



0^A|$) 



A^ ® . . . F N )\*)f M (13) 
= (Fi<g>... F N )® M A\V)® M = 



as a direct consequence of Eqs. ( |TT| ) and (12). 

There are cases in which a state |<£>) can be created 
through local action on a state I*), but where the inverse 
is not possible. In these cases it is meaningful to assert 
that carries entanglement of superior type than the 
the state |$). However, as it is the case for example with 
the GHZ and the W state, there are also different classes 
of entangled states and neither can one state of the former 
class be obtained through local manipulation of a state 
of the latter, nor is the reverse possible. To show such 
distinctness for two states |<J>) and j*) with the help of 



Eq. (12), two operators A\ and A 2 are required with 



= and Ai|$) ^ , but (14) 
A 2 \V) + and A 2 \<$>) = . (15) 

However, if the ranks of all single-subsystem reduced 
density matrices of both states |$) and I*) coincide for 
any subsystem, Eq. (12) is already sufficient for dis- 
tinctiveness: in this case all local Kraus operators can 
be assumed to be invertible, without loss of general- 
ity [5], so that the non-existence of local operators 
(i = l,...,n) with |$) = Fi <g> ... ® F n \&) implies the 
non-existence of local operators Gi (i — 1, n) with 



II. CLASSIFICATION OF ENTANGLED 
MANY-BODY STATES 

Before using the above formalism to identify new struc- 
tures within the set of multi-partite states, let us first 
discuss how prior knowledge can be reproduced, in order 
to gain some insight in this technique. 



Here 'sgn' denotes the signature, i.e. the function that 
yields sgn(n e ) = 1 for any permutation Tl e that can be 
decomposed into an even number of pairwise permuta- 
tions (transpositions), and sgn(n o ) = —1 for any permu- 
tation n o that can be decomposed into an odd number 
of pairwise permutations. Let us inspect the action of A3 
on a three-replica state-vector. If the components of the 
three replicas are pairwise different, then the correspond- 
ing state is mapped on the completely anti-symmetric 
state 
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by A3, i.e. 
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where the signature sgn of a string of numbers is defined 
similarly to the signature of a permutation above, i.e. 
sgn(123) = sgn(231) = sgn(312) = 1 and sgn(213) = 
sgn(321) = sgn(132) = — 1. If i, j and k are not pairwise 
different, then the state-vector is mapped on the zero- 
vector by A3. 

Since three replicas of a state 1*2) contain no term with 
three pairwise different contributions per subsystem, the 
three replicas l^)® 3 are necessarily mapped on the zero- 
vector by the application of A3 on either subsystem. But, 
the three-level entangled state 1*3) behaves differently: 
three replicas contain the component 
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with i =/= j =/= k i, that survive the application of A3, 
such that 



A. The bipartite case revisited 

With Eq. g) we have verified that no entangled bi- 
qubit states ^2} can be created by local action on a 
separable state |*i). To show that in a system of 
higher dimensional constituents a state 1*3) = \/Ai|ll}-|- 
-v/ A2 122) + 1 33) with three finite Schmidt coefficients 
can not be obtained from a state with only two finite 
Schmidt coefficient, i.e. a state with a Schmidt decom- 
position as 1*2), we have to consider three replicas of this 
quantum state. In this framework we can define the pro- 
jector A3 onto the completely antisymmetric (fermionic) 
subspace of a triplicate single-particle Hilbert space, i.e. 
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Together with Eq. (12) this verifies that 1*3) cannot be 
obtained from |^ 2 ) through local manipulation. With ex- 
actly the same argument, one can also show that a state 
\*&n) with N finite Schmidt coefficients can not be ob- 
tained from states with fewer finite Schmidt coefficients, 
simply invoking the projector 
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onto the completely antisymmetric part of nf . Doing 
so, one recovers the full classification of bipartite entan- 
glement in terms of the Schmidt rank |22) . 
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B. GHZ states vs. W states 

The inequivalence of GHZ and W states can be shown 
in a similar manner. Rephrasing the three-qubit tangle 
r [2 3) in terms of permutation operators, one finds that 
r(*) = 16 V / (*| C34 X|*)® 4 with 
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where two consecutive projectors act on the four repli- 
cas of a subsystem as depicted above. The fact that 
A T \GHZ)® A is finite is implied by the finite tangle of the 
GHZ state, but, one also readily convinces oneself that 
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where |£ r ) is the unique eigenvector of A T associated 
with the unit eigenvalue. For this purpose, let us take 
a look at those terms in \GHZ)® A that are not mapped 
on the zero vector by A T . The projection P^ 12 -* p^ 34 ) 
on the first subsystem replica spaces requires the terms 
associated with the first and second replica as well as 
those associated with the third and fourth replica to be 

( 13) (24) 

different. In turn, P_ <g) P_ on the third subsys- 
tem enforces terms associated with the first (second) and 
third (fourth) replica to be different. The only terms in 
\GHZ)® A that satisfy these requirements simultaneously 
are of the form 



(24) 



with i ^ j (i.e. all other terms vanish one by one under 
the application of A T ). What remains to be verified is 
that the term with i = 1, j — 2 does not cancel the term 
with i = 2,j = 1, i.e. 
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Both state-vectors 
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are mapped on the 



singlet state by P_, but with different phases: 



of the first and second subsystem, and in the second and 
fourth replica component of the third subsystem. That 
is, in both of the two cases i = l,j = 2 and i = 2,j = 1, 
the application of A T , which is composed of 6 terms P_, 
yields a sixfold tensor product of \x) with three times 
the negative prefactor — l/v2, and three times the posi- 
tive prefactor l/y/2. That is, in both of these cases, the 
same resulting prefactor —1/8 is obtained, so that these 
two terms add up constructively, i.e. do not cancel each 
other. 

As we discuss in the following, the \W) state behaves 
differently, i.e. A,-|W)® 4 = 0. To verify this, let us intro- 
duce the short hand notations \u>\) = |211), \u>2) = 1 121) 
and = 1 112) . Now, the quadruplication of \W) is 
comprised of 3 4 terms, but we can easily verify that each 
of these terms is mapped on the zero- vector by A T . The 
projector pf 12 - 1 ® p^ 34 ' acting on the replicas of the first 
subsystem, maps all those terms onto the zero- vector that 
do not contain |wi) exactly once in the first two replicas 
and exactly once in both the third and fourth replicas. 
Similarly, the projector P 



(12) 



pi 34 -' , acting on the sec- 



ond subsystem, maps all those terms onto the zero- vector 
that do not contain \u)2) exactly once in the first two 
replica and exactly once in both the third and fourth 
replicas. Thus, the only terms that survive the applica- 
tion of the first- and second-subsystem component of A T 
read 



and 
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However, all these terms arc completely symmetric in the 
third subsystem components, and, therefore, are mapped 
onto the zero- vector by p[ 13) ® P_i 24) on Hf 4 . That is, 
all-together, we have verified that 
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what, together with Eq. (23), veryfies that GHZ states 



cannot be generated through local manipulation of W 
states. Since all reduced density matrices of both the 
GHZ and the W state have full rank, this, in turn, im- 
plies the in-equivalence of these two states. 
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In Eq. (24), both of the terms 
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threefold: 



appear 



appears in the first and second replica 

component of the first and second subsystem, and in the 
first and third replica component of the third subsystem; 

. ) appears in the third and fourth replica component 



C. Beyond qubits 

Having identified the distinctness of GHZ and W 
states, we can now leave the realm of two-level systems 
and identify in a similar fashion inequivalent states in 
higher-dimensional systems. Natural generalizations of 
these states to three-level systems with full rank of all 
reduced density matrices are 
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FIG. 1: On six replicas of a tripartite Hilbert space, with 
three-dimensional subsystems, one can define the operator 
Ast in Eq. (31 1, to show the inequivalence of three-level GHZ 



and W states. It is based on projectors A onto the completely 
anti-symmetric component of triples of local Hilbert spaces. 
Each of the triples that enter Eq. (311 are identified by con- 
necting black lines. 



such a term in the six-fold Wa-state, IW3)® 6 , it neces- 
sarily would need to contain each of the terms |1), |2) 
and 1 3) exactly twice per subsystem. The only terms in 
|W 3 ) that contain the term |3) are |322), |232) and |223). 
If a term in IVF3)® 6 is to contain |3) exactly twice per 
subsystem, i.e. twice in each column in Fig. [l] it must 
solely consist of a six-fold tensor product of the terms 
1 322), 1 232) and 1 223) , each entering twice. Therefore 
any such term in IW3)® 6 does not contain the term |1), 
that is, contains no contribution of the completely anti- 
symmetric three-level state, and therefore is necessarily 
mapped to the zero- vector by A 3r . 

The sixfold \GHZ 3 ) reads |*)® 6 = Ei,i,fc, M> *=i I***}® 
\jjj) ® \kkk) ® \lll) ® \ppp) ® I qqq), but as we discuss in 
the following, only a few of the overall 3 6 terms yield a 
finite contribution after the application of A 3t . Due to 
_4(i23) ac ting on the first subsystem the first three repli- 
cas need to contain three pairwise different components, 
that we label i, j and k. Because of A^ 125 ^ acting on 
the second subsystem, the fifth replica needs to be in the 
state |fc) since the first replicas are in the states \i) and 
\j) . Similarly, the fourth replica needs to be in |j), due 
to the action of A^ 134 ^ on the third subsystem. Finally, 
the sixth replica needs to be in state since .A^ 456 ) 
acts on the first subsystems. Consequently, all terms in 
\GHZ3}® 6 that are not mapped on the zero vector by 
A% T are of the form 



and 



\W*) = 



(|211) + 1121) + 1112) + 1322) + |232) + 1223» . 

(30) 

The operator A T as defined in Eq. ( 22 1 maps the qua- 
druplication of both of these states on something finite, 
i.e. it is not suited to show the inequivalence of these 
two states. To find a suitable operator, we can resort 
to the projector A3 onto the completely anti-symmetric 
component of a triplicate Hilbert space defined above in 
Eq. ( 16 ) in the context of bipartite systems. In terms 
of this projector - that we will simply denote A in the 
remainder - we can define the generalization of A T to 

A 3r = A {123) ® A^ 56) ® A {125) ® A {346) ® A {134) ® *4 (256) 
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(31) 

where the indices ijk of A^ k ' denote the i-th, j-th and 
fc-th replica of the respective single-particle Hilbert space 
as also illustrated in Fig. [TJ 

With this operator, we can explicitly verify that 
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where |£) is the unique eigenvector associated with the 
unit eigenvalue of A 3t . 



As discussed in the context of Eq. (18), only three- 
replica state-vectors with pairwise different components 
survive the application of A. Consequently, if there was 
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What remains to be shown is that the six (3!) different 
terms with pairwise different i, j and k don't add up 
destructively. For this purpose, we have to inspect the 
prefactors as given in Eq. ( 18 1, whichone obtains through 
the application of the six different terms A in A 3r . In 
the first subsystem, there is A acting on the first three 
replicas and on the last three replicas. In these two terms, 
we have 




(34) 



In the second subsystem, A acts on replicas 1, 2 and 5 
and on replicas 3, 4 and 6. In these two terms, we have 
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In the third subsystem, A acts on replicas 1, 3 and 4 and 
on replicas 2, 5 and 6. In these two terms, we have 
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The states [kz), \ks) and \kq) have the same signa- 
ture, whereas the states [m) and \k^} have the opposite 
signature. But, since both signatures appear an even 
number of times, the overall prefactor is always positive, 
so that the different terms of Eq. ( 33 ) add up construc- 
tively after the application of A 3t . With this, we have 
shown the inequivalence of the IW3) and the \GHZ3) 
state. 

However, these two states certainly do not define the 
only classes of entanglement in tripartite three-level sys- 
tems. For example, there is the Aharonov state \x) — 
1/V6((|123) + |231) + |312)) - (|132) + |213) + |321))) 
[23], that also survives A 3t 
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That is, whereas this shows the inequivalence of W3 
and Aharonov states, we need to invoke another sym- 
metry operation to show the inequivalence of GHZ3 and 
Aharonov states. For this purpose, we use the the pro- 
jector 
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onto the completely symmetric (bosonic) part of a trip- 
licate single particle Hilbert space. Similarly to Eq. ( 18 ) 



above, S maps a tri-replica state- vector with components 
|1), 1 2) and |3) on a completely symmetric state 
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With the help of the operators A and S, each acting on 
three replicas of one subsystem, we can now define the 
operator 
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that allows the discrimination of GHZ3 and Aharonov 
states. The triplication of the GHZ 3 state reads 
\GHZ 3 )® 3 = Ef 4 , fe =i l***> ® \333) ® \kkk). All terms 
in this sum with i, j and k not mutually different are 
mapped on the zero-vector by A, since there is no com- 
pletely anti-symmetric three-particle state in less than 
three dimensions. The only terms that are not mapped 
to zero are of the form 



(41) 



with i, j, k pairwise different. The application of A&A&S 
on this terms yields 
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Since all terms enter with a positive prefactor, i.e. add 
constructively, the application of A <& A ® S on the trip- 
licated GHZ state results in something finite: 
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\A) ® \A) ® \S) . (43) 



This is different for the Aharonov state: let us review 
the terms of |x)® 3 m which the terms associated with 
the three replicas of the first and second subsystems si- 
multanuously have pairwise different terms, since only 
those survive the projection A (g> A onto the completely 
anti-symmetric parts. Labelling the replica components 
of the first subsystems \i), \j) and \k) (pairwise different), 
we obtain the two different terms: 
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All factors associated with a single replica component in 
\Xa) have the same signature, and the same also holds 
for all single replica components in \xb)- However, the 
signature of the |x a )-components is different from the 
signature of the |x&)-components, and since the number 
of replicas is odd, \\ a ) enters |x)® 3 with a different sign 
.than \xb)- Similarly, the first subsystem-components of 
\Xa) have the same signature as the second subsystem 
components, and the same holds true for \xb)- That is, 
independently of what this signature is, the contribution 
from the application of A on the first subsystem is can- 
celed by the contribution of A on the second subsystem, 
i.e. 



= ^1-4)® 1-4) • 



(46) 

Since, however, \x a ) and \xb) enter |x)® 3 with a differ- 
ent sign, and the third subsystem components are always 
mapped on the symmetric state \S), these terms cancel 
after the application of A £§> A £3 S, so that 
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Quite surprisingly, we can also see that a biseparable 
state I^bs) = J2i=i V^ilHtp} j with a bipartite entangled 
component of full Schmidt rank, can not be obtained 
from the Aharonov state, since the application of A ® 
A<8)S yields 
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(42) 
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6V6 



|.4)<8> |.4)<g> 



(48) 
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i.e. something finite. If A3 vanishes, i.e. if the state ex- 
hibits only two-qubit entanglement, this component van- 
ishes, and then, indeed, the biseparable state can be ob- 
tained from the Aharonov state. That is, whereas in the 
case of qubits tripartite entanglement is in general su- 
perior to bipartite entanglement, i.e. all bipartite qubit 
entangled states can be obtained from both W and GHZ 
entangled states, this is no longer true for higher dimen- 
sional systems. 



III. OUTLOOK 

The above exemplary cases of incquivalcnt tri-partite 
states of three-level systems and the reproduction of prior 
classification with one and the same recipe, raise confi- 



dence that Eq. (12 1, as originally established in [15] , will 



eventually provide a systematic understanding of the gen- 
eral structure of multi-partite entanglement. Whereas we 
do not rigorously know whether Eq. ( 12 ) is sufficient to 



provide an exhaustive characterization of all incquivalent 
classes of multi-partite entangled states, the following ar- 
guments suggest that this indeed might be the case. On 
the one hand, the rapid growth of the number of different 
permutations with the number M of a state's multiplic- 
ity, which comes in hand with a growing variety of inde- 
pendent terms in Eq. (10 1, results in an enhanced free- 



dom to build suitable operators A. On the other hand, 
any quantitative entanglement property of a quantum 
state, i.e. an entanglement measure, can be expressed in 
form of expectation values of these operators A as well 
[2"T] . Whereas we focussed here exclusively on the classi- 



fication of multi-partite entanglement, our presently de- 
scribed framework bears the potential for much broader 
applications. Consider for example the question of fil- 
tering [25] or entanglement distillation [2B], where local 
Kraus operators Fi that map a mixed state £ to a given 
pure state |^), i.e. F 1 ®. . .®F N T,Fl®. . .<g)F^ = |*)(*|, 
are sought. Here, £ can be the state of an TV-partite 
quantum system, in the case of filtering, or, it can be 
a multiple tensor product of the to-be-distilled state g, 
i.e. £ = g® q , where q subsystems define a single logi- 
cal subsystem as in regular distillation. In this context, 
a relation AZ® M 'A* = 0, but A\V)® M ^ proves the 
non-existence of a filtering- or distillation-protocoll. 

Very naturally, the present framework also facilitates 
the construction of entanglement measures that target 
very specific entanglement properties: considering expec- 
tation values of the symmetry operations A defined here, 
one arrives at scalar functions - just like the tangle |23| . 
that can be defined in terms of Eq. (22), or, general- 
ized concurrences [27] , that can be defined in terms of A. 
Once verified the equality A\&)® A1 — 0, one assures that 
the function (ifr\® M A\Hf)® M is completely insensitive to 
entanglement of the class of |$), what enables the con- 
struction of entanglement measures that access specific 
entanglement properties. 
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